Vortex in axion condensate as a dark matter halo 
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We study vortices in axion condensate on the galactic scale. Such vortices can occur as a result of 
global rotation of the early universe. We study analytical models of vortices and calculate exemplary 
galaxy rotation curves. Depending on the setup it is possible to obtain a vast variety of shapes which 

■ give a good qualitative agreement with observational results. 
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^ '. I. INTRODUCTION 

' The problem of dark matter seems to be a challenge for theoretical and observational investigations. Recent 
\ astronomical observations, like the measurements of cosmic microwave background (CMB) radiation by WMAP 

or measurements of distant supernova type la (SNIa) [3, [Ij, indicate that the Universe, apart from matter, is fulfilled 
^' with a phenomenological fluid with negative pressure (called dark energy) which could be responsible for the current 
1 . acceleration of the universe. While the nature of this energy is still unknown (the cosmological constant A is the most 
O serious candidate), the combination of results from CMB measurements, SNIa data and extragalactic observations 
^ indicate that as much as 2/3 of total energy density of the Universe is in the form of the mysterious dark energy. 
While it dominates the dynamics of the Universe on the large scale, dark matter (matter whose existence has been 
inferred only through its gravity) clearly influences the galactic scale dynamics. Its abundance is given in terms of the 
density parameter ildm = Pdm/Pc where Pc is the critical energy density. For the "concordance" flat ACDM model 
\ we have total matter fim = 0.3; in turn from emission and absorption of photons visible matter is roughly ilvis — 0.04 
and it gives us that dark matter amounts to ildm — 0.26. 

The rotation curves of spiral galaxies give us the strongest evidence for dark matter. This dark matter in its 80% 
is in some form cold and nonbaryonic. Therefore, we can identify only 4% of the total matter content in the present 
universe. If we do not postulate the existence of dark matter it is not possible to explain why at large distance r 
from the centre of a given galaxy, we would find circular velocity ~ GM^^Jr, since visible matter is concentrated 
around its centre. For the observational results of rotation curves see Ref. f4|. However, observations show that Vc 
[ — _ is independent of r at large distances {vc ^ 200km s~^ is its typical value). From numerical simulations of the halo 
. formation we obtain density profiles for both small and large values phaio r^" with a G (1; 1.5) for small and a = 3 
for large distances [1, 0| . 

There are two main candidates for cold dark matter, namely the axion and neutralino 0, @|- this paper we 

■ concentrate on the axion. Axions were originally proposed to solve the strong CP problem in the QCD. If axions have 
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^ . low mass, thus preventing other decay modes, axion theories predict that the universe would be filled with cold Bose- 
Einstein condensate of primordial axions. The axions in this condensate are always nonrelativistic and if the mass 
is about 10~^ eV it would plausibly explain the dark matter problem. There are prospects for direct experimental 
detection of axion. The Axion Dark Matter Experiment (ADMX) ^] searches for weakly interacting axions in the 
dark matter halo of our galaxy. Unfortunately, studies of axion dark matter are not sufficiently sensitive to probe 
the mass regions where axions would be expected. We investigate the possibility of explanation of flat velocity curves 
of spiral galaxies in terms of axion condensate which can be present in the Universe since the Peccei-Quinn phase 
transition. 
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In our investigation we use the Gross-Pitaevski equation in an expanding FRW universe 

where i7(r) is the external potential, is the coupling constant between axions and a{t) is the scale factor. Here we 
describe the condensate by one particle wave function 0(r, t). The circulation in condensate is expressed by 

T = i v-dr = —2711, (2) 
Jc rn 

where I is an integer called topological charge. C denotes any contour around a vortex. When there is no vortex, 
I = and the circulation vanishes. When the condensate is inside a rotating environment the vortex is formed. In the 
interacting condensate vortices with / > 1 are unstable and decay to the vortices with 1 = 1. So in a realistic situation 
we have a net of elementary vortices {I = 1) which are stable. Such nets of vortices are observed in laboratories ^llj . 

We study the possibility that a galactic halo is just such a vortex. We suppose that such a mechanism can be 
realised in the early universe. Namely, in the presence of global rotation, the proposed mechanism yields a huge 
amount of small vortices whose topological charges equal 1. In this paper we consider a singular vortex in the axion 
condensate. 

The question whether rotation is an attribute of the Universe as a whole has been investigated since classical works 
of Lanczos [11], Gamov [ll], Godel 14], Hawking [l^] and recently by Chapline and Mazur [isl. When compared 
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with the CMB anisotropics, the effects of rotation should not be big today '16']. Barrow et al. [l7| showed that the 
cosmic vorticity depends strongly on the cosmological model and for a flat universe the bound on vorticity relative to 
the Hubble parameter at present is uj/H = 2 x 10~^. 

II. FREE AXION CONDENSATE 

First, we must solve equation ([1]) to obtain the density distribution in the vortex. We know that the interaction 
between axions is very week [l9| , so in the first approximation we can omit nonlinear term and external potential in 
equation ([1]). Then we consider, in fact, a free particle equation 



^ 2_9 _ l{l + l) 



{r,t). (3) 



^dt 2a{t) ' ' 2md?(t) 
In case of the de Sitter evolution of universe (a(<) = cxp(iJt)) the solution has the form 

0(f,i) = C • exp + i|ie-2^*|j,(VAr)l^"(0,^), (4) 

where A = 2m/i/?i^. Here jiix) is the spherical Bcssel function related to the ordinary one with 

jiix) ^ ^l^Ji+^ix). (5) 

and Yj" are the standard spherical harmonics. 

Solution ([4]) is non-normalisable. To obtain a normalisable solution we must introduce proper boundary conditions, 
for example (j){r > f) = for some f. But it is not natural in our case, as our Universe is filled homogeneously by 
the axion condensate. Free axions should thus be described approximately by a non-normalisable waves with defined 
energies. Accordingly, our solution describes an axion condensate with defined nonquantized energy. 

We neglected the interaction between axions in the above, but the gravitational interaction on the global scale may 
be important. In the next section we introduce a mean potential to better describe the gravitational interactions of 
axions. 

Our main idea is that small vortices created in the early universe can grow during expansion of the universe. At 
the early stage, the vortex is a quantum object but during expansion it becomes a classical object. We assume that 
during this transition the shape of the vortex is conserved. To calculate how much the vortex h as g rown we need to 
know when its creation took place. To obtain an upper bound for this value we use the relation [l9| 

/ IQi^GeVX 

rua ^ Q^i eV I j 1 , (6) 
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where fa is the energy of Upq{1) symmetry breaking and creation of axions. There is some evidence that ttIq = 10 ^eV 
[20| . for this value of energy is 



/a ~ 6 • lO^GeV. 

Taking the scale factor ao = 1 (at the time of the creation of condensate) we have today 



fa 



6 • lO'^GeV 



3.7K 3.2 • 10-i3GeV 



2 • 10^ 



(7) 



(8) 



So we see that the vortex cannot have grown more than ~ 10^^ times. Now we need to calculate the value of the 
scale factor a inferred from our solutions. Let us examine first the case without the vortex I = 0. The solution is of 
the form 



(f,t) = C(t) 



sin {VXr) 



(9) 



This solution describes a spherically symmetrical halo with the density distribution p{r,t) |(/)(r, f)p. We must 
remember that the coordinate r in not a physical distance, which is R — a-r. For a spherical distribution we calculate 
the velocity rotation curve from the relation 



v{R) = 

where M.{R) is the mass function and is expressed as 



GMiR) 



R 



M{R)^4tt / R''^p{R')dR'. 
Jo 

Figure ([T]) shows the velocity curve for axion condensate with assumed ^/\/a = 0.2 kpc. 
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FIG. 1: The contribution to the galactic velocity curve given by the axion condensate with assumed \f\/a = 0.2 kpc. 

This result describes the observed galaxy velocity curves quite well - we can see the characteristic change to a 
plateau in the rotation curve. For further analysis we choose the point of transition to plateau in the rotation curve 
as our observable, which is estimated to take place for i?trans — 10 kpc {R is the physical value, measured today). 
From our theory we have (for I = 0) 



Rt 



(12) 
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Taking ^ = we have 

Taking (fT^ and (|13p for given i?trans we obtain the relation 

a~(^)-1024. (14) 
When we compare this relation with the value in ([H]) we obtain the restriction for axion velocities 

It is naturally satisfied by axions as gravitationally induced velocities of axions are of the order 

r-) ^ 10"^ (16) 

\ C / grav 

Comparing relations ^ and we can see that in this model, the creation of vortices took place immediately after 
creation of axions. 

Very similar results can be obtained considering the vortex solution with I — \. The velocity curves also exhibit a 
plateau, like the solution with I = 0. The only difference is that we now have the angle dependence p ^ cos^d. This 
produces, like for each vertex solution, distortions from the spherical shape of a dark matter halo. 

Until now we have considered non-normalisable one-particle states. This state can describe the ensemble of axions 
only if they have the same energy. In fact, axions can have general distributions of energy. The dispersion of velocities 
of axions is very small [l9| but it can produce normalised states. Thus, we construct the condensate as a superposition 
of one particle states 

f(r,i)-y dkf{k)Mr,t) = CYr{0,ip)e--^"' J dkf{k)exp(-—e-^"yi{kr). (17) 

It is natural to choose /(fc) as the Gaussian distribution, but in this case we cannot compute the integral explicitly. 
That is why we choose the function 

/(fc) = fc^e-C^A-")', (18) 

which could be considered a shifted Gauss function and leads to analytical expressions. Unfortunately, dispersion and 
shift are now coupled. Assuming this distribution we obtain a normalised wave function 



«'(r,i) 



fc3/2 ^.^r/rof 



2\ TTW'- 



SmHa'^ 



3/4 ' 



(19) 



where 




The velocity curve for this model with arp = 10 kpc is shown in Figure [21 As we can see in this figure, after taking 
a normalisable superposition of states, the velocity curve no longer has a plateau. For higher distances the velocity 
starts to decrease. 

Let us now investigate the length scale of axion distribution. Taking 

— = 1, (21) 
aro 



for R = 10 kpc, we obtain the relation 



1.6 (^)- 10^4, (22) 
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FIG. 2: The contribution to the galactic velocity curve given by the axion condensate with aro = 10 kpc. 



which is equivalent to relation ((TJ]). The case of interest to us is / = 1 because it is a stable vortex configuration. We 
then take 



1 Vg 



■ cos tie 



tip 



and the solution is 



(VAr)2 



(23) 



(24) 



Note, that in the case considered, the density distribution strongly depends on the angle 9, namely p oc cos^ 0. Next, 
we construct a superposition of one-particle states using the distribution function 



f{k) ^ fce-^C^/''^")' 
This leads to the expression for the wave function 



«'(i?,0 = C 



1 1 



erf I -- 
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where 
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a 1 
V2fco ^ 

h 1 



a he fVQ 
\/2) 



exp 



4e?(i + »6) 



cos Oe" 



(25) 
(26) 

(27) 
(28) 



Choosing the realistic values of parameters m — 10^^ eV, H ~ 70 km s^-'^Mpc^^, a = 10^^, uq/c = 10"'^ we obtain 
£_i = 4.5 kpc and ^2 — 3.4 • 10^^'^. Using relation (fTO|) we can now compute the velocity curve. Because relation (fTO|) 
is valid for spherical distributions we compute only the mean velocity curve. This result is shown in Figure [31 



III. GRAVITATIONAL INTERACTION OF AXIONS 



Up to now we have neglected the interaction of axions. But in fact, the dispersion of velocities is a result of 
interaction. This effect was studied in the previous section. Now we include the gravitational interaction by means 
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FIG. 3: The contribution to the galactic velocity curve given by the axion condensate with — 4.5 kpc and ^2 ~ 3.4 -10 . 

of an effective potential. We introduce the average density distribution of axions (p) = {p)o/a'^ which includes both 
condensed and noncondensed parts. Since we work in the Newtonian approximation, we solve the Poisson equation 



aHt) 



(29) 



which in this case leads to the effective potential 



2^2 



U{r) = 



1 muj r 

2 a{t) ' 



(30) 



where 



-7rG(p)c 



(31) 



Now we have to solve the equation 

^fd 3d(t)\ , 1 



2 d + 



Qj,2 J, Qj, 



ir,t) + ^^^mu;\^<P{r,t). (32) 
2 a[t) 



This time it is not possible to separate variables as easily as it was done before, because the right hand side of the 
equation explicitly depends on time. Relative contributions of the right hand side components depend on the factor 
I /a. So we can see that during the evolution, the interaction term becomes more and more important. We cannot 
find the general solution of equation (I32p but we can investigate the boundary solutions. Taking a = 1, we obtain 



(r, a = 1) = — e 
ro 



- ] iFi 

ro 



, 3 (r 

-nr,l + -; — 
2 V^o 



(33) 



where r^^^ = y^muj/h, and i-Fi if is the confluent hypergeometric function of the first kind. The value of the chemical 
potential is quantised as follows 



/i = tiuj ( 2ri,. + ' + 2 



(34) 
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IV. SUMMARY 



We have considered axion condensate as a candidate for dark matter in galactic halos. Condensates in laboratories 
are known to exhibit vortices induced by rotation of the environment. Such vortices are unstable and decay into states 
with the lowest angular momentum. We prolong this picture onto the entire Universe whose global rotation might 
give rise to the rotation of particular galaxies. 

Summing-up, we firstly considered non-normalisable states with I ~ and I = 1 (vortex). It was shown that these 
solutions leads to the flat velocity curves characterised by oscillations in a plateau region. Next, we introduced the 
axion velocity dispersion to obtain well normalised states. Such a velocity dispersion might come from gravitational 
interactions of axions. In this case we calculated contributions to velocity curves which, after reaching a maximum, 
slowly decreases. Finally we mentioned the problem of axion interaction, which requires a further numerical analysis. 

Beside the problem of velocity curves the scenario considered can explain also the galaxy mass problem. As we 
know masses of the spiral galaxies are similar and equal to ~ 10^^ Mq. This can be naturally explained since all 
vortices are similar, what is directly a quantum effect. 
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